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20.  Continued 


Attempts  to  replace  the  Kalman  filter  with  a simple  filter  with  comparable  performance  have 
no.  been  productive.  The  basic  reason  behind  this  difficulty  is  that  accurate  position  and  velocity 
estimates  (obtainable  by  triangulation  from  different  platforms)  require  the  processing  of  position 
and  velocity  covariance  matrices.  Since  both  matrices  must  be  saved  and  updated,  a simple  filter 
does  not  seem  possible. 


TRACKING  FILTERS  FOR  MULTIPLE-PLATFORM  RADAR  INTEGRATION 


INTRODUCTION 

In  the  1960’s,  fleet  exercises  demonstrated  that  many  targets  were  not  detected  by 
radar  operators.  Furthermore  postanalysis  of  video  recordings  of  the  radar  data  revealed 
that  the  radar  return  from  the  targets  was  present  in  the  raw  video.  Some  of  the  reasons 
operators  missed  targets  were  operator  fatigue,  collapsing  of  upper  beams  of  the  3D  radar 
onto  a PPI  display,  jamming,  and  clutter.  Therefore,  to  improve  its  surveillance  capa- 
bility, the  Navy  decided  to  associate  automatic  detection  and  tracking  (ADT)  systems 
with  its  radars.  Specifically  the  SPS-48C  and  RVP  (radar  video  processor)  for  2D  radars 
have  been  approved  for  fleet  operation. 

On  board  most  naval  combat  vessels  there  are  two  kinds  of  surveillance  radars:  2D 
radars  (usually  UHF  band)  and  3D  radars  (S  band).  In  1973  the  Naval  Research  Labora- 
tory [1-4]  and  the  Applied  Physics  Laboratory  (APL)  of  the  Johns  Hopkins  University 
[5]  embarked  on  programs  to  maximize  the  information  aboard  the  vessel  by  integrating 
the  radar  data  from  different  radars  into  a single  track  file.  The  benefits  of  such  a system 
would  be  increased  track  life  and  on-line  redundancy.  At  present  the  SYS-l-D  system,  an 
operational  automatic  detection  and  integrated  tracking  (ADIT)  system  developed  by 
APL,  is  scheduled  to  be  tested  at  sea  in  early  1977. 

As  an  extension  of  the  concept  of  integration  more  information  can  be  obtained  by 
combining  .he  radar  data  from  several  platforms  into  a system  track  file.  In  addition  to 
the  advantages  of  an  ADIT  system,  a platform-to-platform  integration  system  will  have 
inherent  antijamming  protection  because  of  its  frequency  arid  spatial  diversity.  Some  of 
the  basic  ideas  and  problems  associated  with  such  a system  are  discussed  in  NRL 
Memorandum  Report  3404  [6].  In  the  current  report  some  tracking  filters  that  can  be 
used  for  platform-to-platform  tracking  are  discussed. 

The  basic  tracking  philosophy  and  two  tracking  algorithms  (Kalman  filter  and  an 
uncoupled  maximum-likelihood  filter)  are  discussed  in  the  next  section.  The  accuracy  of 
the  two  filters  is  compared  by  using  a Monte  Carlo  simulation.  The  simulation  is  discussed 
in  the  third  section,  and  the  results  are  given  in  the  fourth  section.  The  conclusions  are 
made  in  the  fifth  section. 


TRACKING  FILTERS 

Since  the  raw  detections  contain  the  maximum  amount  of  information,  communica- 
tion bandwidth  restrictions  will  be  ignored,  and  it  will  be  assumed  that  all  detections  will 
be  used  to  update  a track.  All  detections  do  not  contain  the  same  amount  of  informa- 
tion [6] : one  radar  could  be  more  accurate  than  another,  or  a second  detection  could 
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immediately  follow  (in  time)  a previous  detection.  Consequently  what  is  required  is  a 
method  of  filtering  the  data  (before  transmission  over  the  communication  channel)  to 
obtain  the  detections  containing  the  most  information.  Then  only  these  detections 
would  be  used  to  update  the  track.  However,  since  the  question  of  filtering  the  detec- 
tions is  still  open,  it  is  assumed  in  this  report  that  all  detections  are  used  to  update  the 
track. 

The  first  filter  considered  is  the  Kalman  filter  in  xy  coordinates.  Since  the  storage 
and  computation  requirements  are  somewhat  large,  an  uncoupled  maximum -likelihood 
filter  is  also  considered. 


Kalman  Filter 

The  Kalman  filter  is  a recursive  filter  which  minimizes  the  least-square  error.  The 
state  equation  in  xy  coordinates  [7],  which  in  our  case  represents  the  equation  of 
motion,  is 


X(t  + l)  = *(t)X(t)  + r(tM(t),  (1) 


~x(t)~ 

"l  T 0 0~ 

— T2  0 
2 

x(t) 

0 10  0 

T 0 

~ax(f) 
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, ‘t’(f)  = 
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with  X{t)  being  the  state  vector  at  time  t consisting  of  position  and  velocity  components 
x(t),  x{t),  y(t),  and  y(f),  t + 1 being  the  next  observation  time,  T being  the  time  between 
observations,  and  ax(t)  and  (t)  being  random  accelerations  whose  covariance  matrix  is 
Q(t).  The  observation  equation  is 

Y(t)  = M(t)X(t)  + V(t),  (2) 


where  Y(t) 


a Xm 
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mt)  = 


10  0 0 
0 0 10 
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with  Y(t)  being  the  measurement  at  time  t consisting  of  positions  xm(t)  and  ym  (f)  and 
V(t ) being  zero  mean  noise  whose  covariance  matrix  is  R(t). 

The  problem  is  solved  recursively  by  first  assuming  the  problem  is  solved  at  time 
t - 1.  Specifically  it  is  assumed  that  the  best  estimate  X(t  - l|t  - 1)  at  time  t - 1 and 
its  error  covariance  matrix  P{t  - lit  - 1)  are  known,  where  the  circumflex  signifies  an 
estimate  and  X(f|s)  signifies  that  X(t)  is  being  estimated  with  observations  up  to  F(s). 
The  six  steps  involved  in  the  recursive  algorithm  are  as  follows: 
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Step  1.  Calculate  one  step  prediction, 

X(t|t-l)  = 4>(t-l)X(f-l|t-l); 

Step  2.  Calculate  the  covariance  matrix  for  one  step  prediction, 

P(f|f  - 1)  = 4>(t  - l)P(t  - lit  - l)$(f  - 1)  + T(t  - 1)  Q(t  - 1)  T(t  - 1); 
Step  3.  Calculate  the  prediction  observation. 


y(f|t-l)=M(0*(t|t-l); 


Step  4.  Calculate  the  filter  gain, 


A(t)  = P(f  It  - l)Af(t)  [Af(t)P(tU  - 1)  M(t)  + R(t)]  ■ 


Step  5.  Calculate  a new  smoothed  estimate, 


X(t|t)  = X(t|t  - 1)  + A(f)  [y(t)  - 7(t|t  - 1)] ; (7) 

Step  6.  Calculate  a new  covariance  matrix, 

P(tlt)  = (/-A(t)M(t)]P(t|t-l).  (8) 

In  summary,  starting  with  an  estimate  X(t  — 1 It  — 1)  and  its  covariance  matrix  P(f  — 1 |t  - 1), 
after  receiving  a new  observation  Y(t)  and  calculating  the  six  quantities  in  the  recursive 
algorithm,  a new  estimate  X(t|t)  and  its  covariance  matrix  P(t|f)  are  obtained. 

For  the  Kalman  filter  in  xy  coordinates,  the  measurement  covariance  matrix  R(t)  is 
a function  of  the  radar-target  geometry.  Letting  (at  time  t)  rt  and  6t  be  the  range  and 
azimuth  of  the  target  with  respect  to  the  radar  (with  the  azimuth  angle  being  measured 
counterclockwise  from  the  x axis),  the  elements  of  the  covariance  matrix 


*(0  = 


°2y(t) 


o^(t)  = of  cos2  6t  + rfo%  sin2  6t, 

(Jy(t)  = o2  sin2  Of  + rfog  cos2  6t, 
and 

o2v(t)  = [a2  - rfog]  sin  0t  cos  9t, 

where  a2  and  a2  are  the  variances  of  the  range  and  azimuth  measurement  errors 
respectively. 
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Fig.  1 — Typical  contours  of  equal  probability  for  the 
position  of  a target  that  is  detected 


For  the  assumption  that  the  measurement  errors  are  Gaussian,  typical  contours  of 
equal  probability  are  shown  in  Fig.  1.  If  or  = 100  m,  ae  = 0.3°,  r = 2 X 105  m,  and 
6 = 90°,  the  corresponding  covariance  matrix  is 


R 


106  0 " 
0 104 


(13) 


Thus  the  radar  measures  the  y coordinate,  corresponding  to  range,  rather  accurately  and 
the  x coordinate,  corresponding  to  azimuth,  rather  inaccurately.  If  a Kalman  filter  is 
used  for  the  radar  geometry  shown  in  Fig.  1,  the  covariance  matrix  P(t|f)  will  be  of  the 
form 


P(flt)  = 


'H(t) 

0 


0 " 
L(t) 


(14) 


where  H(t)  and  L(f)  are  two-by-two  matrices  and  the  terms  in  H are  approximately  100 
times  greater  than  the  corresponding  terms  in  L.  The  filter  gain  A (t)  is  of  the  form 


Mt)  = 


ax(t)  0 

MO  o 

0 ay(t) 

0 /3y(t) 


(15) 
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MEASUREMENT  1 


A 


MEASUREMENT  2 


\J 


• RADAR  2 


• RADAR  1 

Fig.  2 — Geometry  of  two  radars  and  of  where  they  detect  a target 


and  furthermore  ax(f)  o ’y(t)  and  Px(t)  **  Py(t).  The  filter  gains  would  be  identical  if 

the  target’s  range  change  were  only  a fraction  of  its  initial  range  (the  x -measurement 
variance  r2a2  must  remain  relatively  constant). 

If  at  some  instant  in  time  a second  radar  with  a different  aspect  angle  detects  the 
target,  the  situation  changes  significantly.  For  example,  if  as  shown  in  Fig.  2 a second 
radar  detects  the  target  at  an  azimuth  of  180°  at  the  same  range  with  the  same  accuracy 
as  the  first  radar,  its  measurement  covariance  is  of  the  form 


*(*  + !)  = 


104 

0 


0 

106_ 


(16) 


The  filter  gain  A(t  + 1)  is  of  the  form 


A(f  + 1)  = 


”£**(' + D 
Px(t  + 1) 
0 
0 


0 ~ 
0 

«y(*+l) 

Py(t  + 1)_ 


(17) 


and  ax(t  + 1)  58  1,  Px(t  + 1)  > Px(t),  o 'y(t  + 1)  0,  and  Py(t  + 1)  « 0.  That  is,  initially 

radar  1 was  providing  an  accurate  y measurement  and  an  inaccurate  jc  measurement. 

When  radar  2 provided  an  accurate  x measurement  and  an  inaccurate  y measurement,  the 
accurate  measurement  is  given  a weighting  of  1,  and  the  inaccurate  measurement  is  given 
a weighting  of  0. 


5 


Mn&HRRBBRMMMMOMi 


TRUNK  AND  WILSON 


As  a simple  numerical  example  assume  that  the  track  is  initiated  by  two  detections, 
4 seconds  apart,  made  by  radar  1 in  Fig.  2.  Four  seconds  later  radar  1 detects  the  target 
and  makes  the  first  update  of  the  track  with  the  Kalman  filter  gains 


4(1)  = 


5/6 

0 

1/8 

0 

0 

5/6 

0 

1/8 

(18) 


If  again  radar  1 detects  the  target  4 seconds  later,  the  gains  are  reduced  to 


4(2)  = 


7/10 

0 

3/40 

0 

0 

7/10 

0 

3/40 

However,  if  radar  2 makes  the  next  detection  4 seconds  later,  the  gains  are 

1 0 


4(3)  = 


1/12  0 
o nr2 
o nr3 


(19) 


(20) 


The  gain  for  the  x velocity  rose  only  to  1/12  because  the  old  x measurement,  which  is 
used  to  calculate  the  x velocity,  is  rather  inaccurate.  However,  when  a second  detection 
from  radar  2 is  made,  Px(")  will  rise  considerably,  since  an  accurate  x velocity  can  now 
be  estimated. 


Kalman  Filter  with  a Turn  Detector 

The  Kalman  filter  is  the  optimum  filter  as  long  as  the  target  trajectory  obeys  the 
state  equation  (1),  which  describes  a straight-line  trajectory  with  random  perturbations 
(the  random  perturbations  bound  the  filter  gains  away  from  zero).  However,  when  the 
target  maneuvers,  the  maneuver  must  be  detected  and  the  error  covariance  matrix  must 
be  increased.  In  this  study  the  error  criterion  is 

£ = [X(t|t  - 1)  - X(t|t)]  M [M  P(t\t  - DM]"1  M [X(tlt  - 1)  - X(t|t)j 

+ [7(f)  -AfX(tlf)]  R(t)-1  [7(f)  -MX(t|t)].  (21) 
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This  error  is  the  squared  Mahalanobis  distance  from  the  smooth  position  MX{t\t)  to  the 
predicted  position  MX(t\t  - 1)  plus  the  squared  Mahalanobis  distance  from  the  smooth 
position  MX(t\t)  to  the  measured  position  7(f).  The  Mahalanobis  distance  differs  from 
the  Euclidean  distance  by  using  a covariance-matrix  kernel  instead  of  an  identify  matrix. 

When  the  error  E is  greater  than  a threshold  (which  in  this  study  was  set  to  E = 16, 
corresponding  for  example  to  covariance  matrices  that  are  diagonal  and  smooth  coordi- 
nates that  differ  from  the  predicted  and  measured  positions  by  twice  the  standard 
deviation),  the  error  covariance  matrix  P(t  - Ilf  - 1)  is  increased  and  a new  smooth 
position  MX(t\t)  is  calculated.  Increasing  P(t  - l|t  - 1)  causes  the  new  position^ estimate 
MX(t\t)  to  be  closer  to  the  measurement  Y(t)  and  further  from  the  prediction  X(t|t  - 1). 
Since  P(t\t  - 1)  increases  when  P(t  - l|t  - 1)  increased,  this  increase  in  P(t  - 1 It  - 1)  will 
always  cause  E to  decrease.  This  procedure  is  repeated  until  E is  less  than  the  threshold. 
Specifically  terms  Pn,  P13,  P31,  and  P33  are  increased  by  s/T ; terms  P12,  P14,  P2i,  P23, 
P32,  P34,  P40,  and  P43  are  increased  by  F\  and  terms  P22,  P24,  P42,  and  P44  are 
increased  by  F2.  (In  this  study  F = 1.5”,  where  n is  the  number  of  consecutive  covariance 
matrix  increases.)  The  position  covariance  elements  are  not  increased  as  much  as  the 
velocity  elements  because  of  coupling;  that  is,  an  uncertainty  in  predicted  position  is  due 
not  only  to  the  uncertainty  in  the  last  position  but  also  in  the  velocity.  In  a real  system 
the  track  should  also  be  bifurcated  when  a large  error  is  encountered. 


Computational  Requirements  for  the  Kalman  Filter 

Since  the  computational  load  in  performing  the  Kalman  filter  (equations  (3)  through 
(8))  appears  formidable,  care  was  taken  to  minimize  the  number  of  calculations  (additions 
and  multiplications).  The  Kalman  filter  requires  50  additions  and  61  multiplications. 
These  numbers  take  into  account  the  symmetry  of  covariance  matrices  and  the  simplicity 
of  certain  matrices;  for  example,  HP(t\t  - l)ff  removes  four  elements  from  a tour-by-four 
matrix  and  consequently  requires  no  additions  or  multiplications. 

The  storage  requirement  of  the  Kalman  filter,  in  addition  to  the  positions  and 
velocities,  is  the  ten  unique  elements  of  the  covariance  matrix.  In  an  effort  to  reduce 
the  computational  load  and  the  storage  requirements,  a modified  maximum-likelihood 
approach  was  used. 


Modified  Maximum-Likelihood  Filter 

In  this  subsection,  it  will  be  assumed  that  the  predicted  and  measured  variables  are 
independent  and  Gaussian.  If  the  position  and  velocity  variables  were  considered  jointly, 
the  maximum-likelihood  method  would  yield  the  Kalman  filter.  Consequently,  to  reduce 
the  complexity,  position  and  velocity  are  considered  separately. 


The  joint  density  of  the  predicted  position  Xp  and  measured  position  Xm  is 


P(Xp,Xm)  = 


-1, 


r(\l2)(xp-»)Kp(Xp-n ) 


yi2 wrp 


-(l/2)(*m-M)Km  (*m-v ) 


s/WK„ 


, (22) 


I 
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where  Kp  and  Km  are  the  predicted-position  and  measured-position  covariance  matrices 
respectively.  The  maximum-likelihood  estimate  of  the  position  p is  that  value  of  p which 
maximizes  (22).  Taking  the  partial  derivative  of  the  log  of  (22)  with  respect  to  p yields 
the  expression 

ICpUXp  -p)  + K^(Xm  -p).  (23) 

When  (23)  is  set  equal  to  0 and  solved  for  p.  the  maximum-likelihood  estimate  is  obtained 
as 

p * <Kp'  ♦ Kj  r»  (Kp'.Yp  ♦ Km  Xm  ).  (24) 

It  is  straightforward  to  show  that  the  covariance  of  p is 

Ku  *E(p£  ♦/Cp1)'1.  (25) 

Using  the  new  position  estimate  p and  the  old  smoothed  position,  a new  velocity 
estimate  can  be  obtained.  Then  the  new  velocity  estimate  can  be  combined  with  old 
velocity  estimate  using  an  equation  with  the  identical  form  of  (24). 

In  summary,  the  problem  is  originally  solved  for  a given  instant  in  time.  Specifically 
the  smoothed  position  X,  and  smoothed  velocity  Vt  and  their  corresponding  covariance 
matrices  Kx  and  /CJ1  (inverse)  are  available  After  T seconds  a new  position  Xm  is 
measured;  and  the  new  estimates  and  covariam  s are  found  by  calculating  the  eight 
quantities  in  the  following  algorithm: 

Step  1.  Calculate  the  predicted  position  Xp, 

X p = X,  * V t T.  (26) 

Step  2.  Calculate  the  covariance  matrix  (assuming  .Yf  and  Vs  are  independent)  for 
prediction, 


Kp  = Kx  * T2  Kv;  (27) 

Step  3.  Calculate  the  covariance  matrix  for  the  smoothed  position  estimate, 

« (Kp1  ♦ r1 ; (28) 

Step  4.  Calculate  the  smoothed  position  estimate. 

p = KM  (Xp1  Xp  Xm)\  (29) 

Step  5.  Calculate  the  new  velocity  estimate, 

VN=Cp-X,)IT ; (30) 
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Step  6.  Calculate  the  covariance  matrix  (assuming  p and  Xs  are  independent)  for  the  new 
velocity  estimate, 


Kn  = (KM  + Kx)/T2; 


(31) 


Step  7.  Calculate  the  inverse  covariance  matrix  for  the  smoothed  velocity  estimate, 

K:1  =(Kv1 +K^y,  (32) 

Step  8.  Calculate  the  smoothed  velocity  estimate, 

V = Ki(K;1  Vs+Ktf  VN).  (33) 


These  eight  steps  complete  the  cycle:  (29)  is  the  smoothed  position  estimate,  (28)  is  its 
covariance  matrix,  (33)  is  the  smoothed  velocity  estimate,  and  (32)  is  its  covariance 
matrix. 


Modified  Maximum-Likelihood  Filter  with  a Turn  Detector 


The  turn  detector  is  the  same  one  which  was  used  for  the  Kalman  filter.  In  terms  of 
the  new  parameters  the  error  criterion  is 


E = (Xp  - p)  Kp1  (Xp  - p)  + (Xm  - p)  tf'1  (Xm  - p). 


(34) 


When  the  error  exceeds  16,  the  covariance  matrices  Kf  and  Kv  are  both  increased  by  the 
factor  F (presently  F = 1.25).  The  covariance  matrices  are  continually  increased  by  F 
until  E < 16.  Although  there  is  no  reason  for  increasing  Ks  and  Kv  by  the  same  factor, 
this  procedure  has  led  to  good  results. 


Computational  Requirements  For  the  Modified 
Maximum-Likelihood  Filter 


Taking  into  account  the  symmetry  of  the  covariance  matrices,  the  maximum-likelihood 
method  (equations  (26)  through  (33))  requires  37  additions  and  59  multiplications.  The 
storage  requirements,  in  addition  to  the  positions  and  velocities,  are  the  six  unique 
covariance  elements  of  the  two  covariance  matrices  Ks  and  Kv.  Comparing  these  results 
with  those  of  the  Kalman  filter,  it  is  seen  that  the  maximum-likelihood  method  provides 
little  computational  advantage.  This  is  because  the  calculations  are  essentially  repeated 
twice  for  the  maximum-likelihood  method:  once  for  the  smooth  position  and  once  for 
the  smooth  velocity.  All  attempts  to  simplify  the  filter  by  using  only  one  covariance 
matrix  have  failed.  The  reason  is  demonstrated  by  the  example’s  gains  given  in  (18),  (19), 
and  (20).  In  the  example,  when  the  system  receives  its  first  measurement  from  radar  2, 
the  position  gain  is  set  close  to  1.  However  it  is  not  until  the  second  measurement  from 
radar  2 that  the  velocity  gain  will  rise  significantly.  Since  there  may  be  a considerable 
time  delay  in  obtaining  this  second  measurement,  two  covariance  matrices  must  be  saved. 
The  need  to  save  two  covariance  matrices,  with  the  cross  terms  supplying  the  crucial 
information  for  triangulation  of  a target,  has  aborted  all  attempts  to  simplify  the  tracking 
filter. 
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Radar  Geometry 

The  simulation  involves  two  platforms.  The  first  ship  is  centered  at  the  origin  of 
an  xy  coordinate  system  (xx  = yt  =0)  whose  x axis  is  through  the  bow  of  the  ship.  The 
second  ship  is  at  coordinates  (x2,  >2)  and  *s  oriented  parallel  to  the  first  ship.  Azimuth 
is  measured  as  usual  in  an  xy  coordinate  system:  target  at  (x,  y)  with  respect  to  radar  i 
has  azimuth  angle  0,  = tan-1  [(y  - y,)/(x  - x,)].  The  rolling  and  pitching  of  the  ships  are 
assumed  to  be  periodic  and  are  given  by 


R,(t)  = Rm  sin  + 7f ) 


(35) 


and 


P,(t)  = PM  sin 


( 2irt  „ \ 

* +£l)' 


(36) 


where  Rm  and  Pm  are  the  maximum  roll  and  pitch  angles,  Tr  and  T/>  are  the  roll  and 
pitch  periods,  which  are  independent  and  uniformly  distributed  between  10  and  12 
seconds,  and  7,  and  £,•  are  independent  phase  angles  uniformly  distributed  on  0 to  2ir. 


Target  Trajectory 


The  target  flight  profile  is  specified  by  an  initial  range,  altitude,  azimuth,  speed,  and 
heading.  The  target  proceeds  along  this  heading  until  time  ts,  when  the  target  starts  a 
counterclockwise  turn  at  a specified  g value.  At  te  the  target  ends  the  turn  and  proceeds 
along  a straight  line  at  its  present  heading.  The  elevation  angle  of  the  target  is  calculated 
by  letting  x,  y,  and  z be  the  target  coordinates  on  a flat  earth  and  x,-,  y,-,  and  hr  be  the 
radar  coordinates.  Then,  if  the  4/3  radius  of  the  earth  is  denoted  by  Re  and  if  the 
notation 


and 


/ ( hr\ 

■{°-5rJ  -Gr) 


)2jl/2 

(37) 

.1/2 

(38) 

is  introduced,  then  the  elevation  angle  e of  the  target  is  given  [8]  by 


t>  = tan 


-1 


(».-) 


(39) 
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Target  Detection  and  Estimation 

The  scan  time  for  radar  1 is  uniformly  distributed  between  3.9  and  4.1  seconds,  and 
the  scan  time  for  radar  2 is  uniformly  distributed  between  4.9  and  5.1  seconds.  The 
initial  time  radar  1 passes  over  the  target  is  t = 0,  and  the  initial  time  radar  2 sweeps  past 
the  target  is  uniformly  distributed  between  zero  and  the  scan  time  of  radar  2. 

The  question  of  whether  or  not  a target  is  detected  is  resolved  by  first  calculating 
the  target  signal-to-noise  ratio  (S/N).  Using  Blake’s  model  [9] , the  multipath  propagation 
factor  F is  calculated.  Then  the  S/N  is  calculated  using 


(40) 


S/N  = 20  oT 


where  aT  is  the  target  cross  section  in  square  meters,  RD  is  the  range  where  the  probabil- 
ity of  detection  PD  is  0.9  and  the  probability  of  false  alarm  Pfa  is  10-6  for  a 1-square- 
meter  target,  and  R is  the  target  range.  Using  curves  from  Robertson  [10],  one  obtains 
the  PD  at  = 10'6  for  the  calculated  S/N.  A uniformly  distributed  random  number 
between  0 and  1 is  generated,  and  if  the  number  is  below  PD.  the  target  is  declared  to  be 
detected. 


If  the  antenna  is  unstabilized  and  if  a,  and  e,  are  the  true  azimuth  and  elevation 
angles  of  the  target,  the  measured  angles  are  [11] 


am  <0  = tan' 


|" sin  a,  cos  P,  + (cos  a,  sin  P,  + tan  ei  cos  P,)  sin  P,  1 


cos  o,  cos  P,  - tan  ei  sin  P,- 


+ e. 


(41) 


and 


em  (i)  = sin-1  [cos  e,  cos  a,  sin  P,  + sin  e,  cos  Pf)  cos  P,-  - cos  e,  sin  a,  sin  P,]  + ei2,  (42) 

where  e,!  and  e|2  are  independent  Gaussian  random  variables  with  variances  a2  and  a2 
respectively.  Letting  U and  V be  independent  uniformly  distributed  random  variables 
between  0 and  1,  Gaussian  variables  with  mean  zero  and  variance  a2  can  be  generated  by 

e = a(-  2 log  U)1/2  cos  (27rV).  (43) 

If  one  measured  only  am(i)  and  not  em(i),  one  would  have  large  azimuth  errors.  For 
instance,  if  Rm  = 10°,  Pm  = 5°,  and  e = 15°,  the  azimuth  error  can  be  as  large  as  5°  even 
though  oa  - 0.5°  [12].  However,  if  em(i)  is  measured  and  P,  and  P,  are  known,  the 
measurements  am(i)  and  em(i)  which  are  relative  to  the  deck  plane  of  the  ship  can  be 
rotated  into  a system  whose  jcy  plane  is  the  plane  of  the  ocean.  These  equations,  which 
were  derived  by  George  as  cited  in  Ref.  13,  are 


a'm(i)  = tan-1 


- sin  P,  sin  em  ( i ) + cos  P,  sin  am  (i)  cos  em  (i) 
cos  P,  cos  am  (i)  cos  em  (/)  + IV  sin  P, 


(44) 
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e'm  (i)  = sin-1  [sin  P,  cos  am  (i)  cos  em  (i)  + W cos  P,] , (45) 

where 

W = cos  Ri  sin  em  ( i ) + sin  P,  sin  am  (i)  cos  em  (j). 

Radars  involved  in  platform-to-platform  radar  integration  need  elevation  information  [6] 
to  transmit  useful  information  from  one  platform  to  another.  Thus,  elevation  angles  are 
available  to  perform  the  appropriate  corrections. 


MONTE  CARLO  RESULTS 

In  this  section  results  of  the  Kalman  and  maximum-likelihood  filters  are  obtained 
for  the  radar  geometry  and  target  trajectories  shown  in  Fig.  3.  The  radar  coordinates  are 
(0,  0)  and  (60,  -60)  km  and  the  radar  heights  jure  23  m.  The  target-trajectory  parameters 
are  given  in  Table  1. 

It  will  be  assumed  that  both  radars  have  the  same  accuracies.  Specifically  the 
standard  deviations  of  the  range,  azimuth,  and  elevation  measurements  are  or  = 150  m, 
oa  = 0.3°,  and  oe  = 0.3°. 

The  performance  measure  for  the  tracking  filter  will  be  a modified  RMS  velocity 
error  defined  by 


^RMS 


Vx(i)]2  + [Vy(i)  - Vy(i)} 


(46) 


Fig.  3 — Geometry  of  two  radars  involved  in  platform-to-platform 
radar  integration  and  of  three  cases  of  target  trajectories 
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Table  1 — Target  Trajectory  Parameters  (shown  in  Fig.  3) 


Case 

Range 

(n.mi.) 

Azimuth 

(deg) 

Heading 

(deg) 

Speed 

(m/s) 

Altitude 

(m) 

Cross 

Section 

(m2) 

Target  Maneuver 

— 

Roll 

Rm 

(deg) 

Pitch 

PM 

(deg) 

Start 

(s) 

End 

(s) 

Turn 

(g’s) 

A 

100 

45 

-90 

300 

15000 

0.5 

300 

300 

3 

15 

5 

B 

120 

45 

-135 

300 

4500 

0.5 

100 

111.0 

3 

15 

5 

C 

100 

0 

0 

300 

4500 

1.0 

100 

132.7 

3 

15 

5 

where  Vx(i)  and  Vy ( i)  are  the  true  velocities,  Vx(i)  and  Vy(i)  are  the  estimated  velocities, 
and  the  asterisk  indicates  that  the  sum  excludes  points  between  ts  and  te  + 10.  The 
samples  at  the  turn  are  excluded  because  during  a turn  the  errors  are  large  and  would 
dominate  VRMS . During  a turn  what  is  important  is  the  ability  to  detect  the  turn.  Since 
both  filters  have  the  same  turn  detector  (comparison  of  (21)  to  16),  the  filters  will  have 
the  same  performance  in  this  area.  In  the  simulation  the  radar  platforms  are  assumed  to 
be  gridlocked. 

The  simulation  was  run  five  times  for  the  two  filters  and  three  target  trajectories, 
and  the  average  VRMS  is  given  in  Table  2.  As  one  would  expect,  the  Kalman  filter  pro- 
vides the  more  accurate  track.  Comparison  of  case-by-case  results  show  that  the  Kalman 
filter  reacts  more  rapidly  when  a poor  initial  estimate  is  made.  For  example,  the  velocity 
estimates  for  the  third  run  of  case  B are  shown  in  Table  3.  The  initial  velocities  are  the 
same,  the  initialization  algorithm  being  the  same.  As  can  be  seen  from  Table  3,  while 
the  velocities  from  both  filters  are  approaching  the  true  values  (V*  = Vy  = -215  m/s),  the 
velocity  from  the  Kalman  filter  is  approaching  more  rapidly. 

Probably  adjustments  could  be  made  to  quicken  the  convergence  of  the  maximum- 
likelihood  algorithm.  However,  since  the  maximum-likelihood  filter  is  almost  as  compli- 
cated as  the  Kalman  filter,  no  changes  were  made,  because  the  Kalman  filter  would  be 
used  rather  than  the  maximum-likelihood  filter. 


CONCLUSIONS 

The  Kalman  filter  is  the  optimum  tracking  filter  (with  respect  to  the  mean-square 
error)  regardless  of  whether  or  not  the  radar  detections  are  made  from  single  or  multiple 
platforms.  The  performance  (specifically  the  RMS  velocity  error)  of  the  Kalman  filter 
for  two  platforms,  separated  by  46  n.mi.,  has  been  calculated  for  various  target  trajectories. 
An  error  criterion  involving  the  squared  Mahalanobis  distances  between  the  smooth  and 
predicted  positions  and  the  smooth  and  measured  positions  is  used  to  detect  target 
maneuv  ts.  After  a turn  has  been  detected,  the  covariance  matrix  is  increased  (in  steps) 
until  the  error  criterion  is  below  a critical  value. 

Attempts  to  find  a simple  filter,  with  good  performance,  have  not  been  productive. 
The  maximum-likelihood  filter  obtained  by  arbitrarily  decoupling  the  position  and 
velocity  estimates,  although  obtaining  acceptable  performance,  is  almost  as  complicated  as 
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Table  2 — Average  VRMS  (Eq.  (46)) 
For  Five  Runs  of  Each  of  the  Cases 
in  Fig.  3 


Case 

Average  VRMS  (m/s) 

Kalman 

Filter 

Maximum- 

Likelihood 

Filter 

A 

16 

23 

B 

32 

60 

C 

22 

32 

Table  3 — Velocity  Estimates  During  the 
Third  Run  of  Case  B 


Velocity  Estimates  (m/s) 

Maximum- 

Time 

Kalman  Filter 

Likelihood 

Filter 

warn 

| 

| 

D 

12.1 

23 

-451 

23 

-451 

16.2 

-125 

-310 

-30 

20.2 

-154 

-274 

-57 

24.3 

-174 

-248 

-82 

-341 

27.4 

-162 

-261 

-105 

-316 

28.3 

-161 

-264 

-104 

-317 

32.4 

-164 

-258 

-118 

-301 

32.4 

-170 

-251 

-118 

-301 

36.4 

-189 

-232 

-146 

-272 

the  Kalman  filter  and  consequently  would  not  be  used.  The  basic  reason  behind  this 
difficulty  is  that  accurate  position  and  velocity  estimates  (obtained  by  triangulation) 
require  the  processing  of  position  and  velocity  covariance  matrices.  Since  both  matrices 
must  be  saved  and  updated,  a simple  filter  does  not  seem  possible. 

In  summary,  the  Kalman  filter  with  a turn  detector  should  be  used  as  the  tracking 
filter  for  radar  detections  from  multiple  platforms. 
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